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I. INTRODUCTION 


Photon-neutrino processes are potentially of interest in astrophysics and cosmology. 
Four-leg elastic scattering processes are suppressed by powers of uj/Miy, where oj is the 
centre-of-mass energy of the collision and Mw the W boson mass [^. Their cross sections 
are therefore too small to be of astrophysical interest |^. On the other hand, five-leg 
processes involving two neutrinos and three photons, such as 

yz/ —>■ 77Z/ 

77 —>■ yz/P 

z/P ^ 777 , ( 1 ) 

are not negligible. In fact, the extra a in the cross section is compensated by an interchange 
of the uj/Mw suppression by an a;/me enhancement |^. 

In ref. [Q , Dicus and Repko derived an effective Lagrangian for the above hve-leg photon- 
neutrino interactions, by substituting one photon with a neutrino current in the Euler- 
Heisenberg Lagrangian that describes the photon-photon scattering |^. Some of the for¬ 
mulae reported in ref. |Q were in disagreement with the results derived earlier by Van Hieu 
and Shabalin [^]. To settle this question, in a recent work |0] we computed the hrst and the 
second process in eq. (i, in the framework of the effective theory, conhrming the results 
reported in ref. |Q. We also justihed the derivation of the hve-leg effective vertex starting 
from the Euler-Heisenberg Lagrangian. 

The effective approach gives reliable results for energies below the threshold for e''“e“ 
pair production, while its extrapolation to energies above 1 MeV, interesting to study, for 
example supernova dynamics, is suspect. Therefore, an exact calculation of the processes in 
eq. (|l|) is important in order to dehnitively assess their role in astrophysics and the range of 
validity of the effective theory. Such a calculation, assuming massless neutrinos, is the main 
ingredient of this paper. The range of energy in which the above reactions are relevant is 
well below the W mass, so that we treated the neutrino-electron coupling as a four-Fermi 
interaction. 

Very recently, parallel work in the same direction has been carried out by Dicus, Kao and 
Repko 0] , so that we a chance to compare our numerical results, finding complete agreement 
between the two independent calculations. 

An additional problem is searching for simple approximations to the exact result. We 
tried the following two approaches: 

• we expanded the exact amplitude in powers of cn/mg, to derive the first correction 
term to the effective theory; 
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• we fitted the curves obtained with the complete calculation. 

It turned out that, as expected in ref. 0, an expansion in inverse powers of the electron 
mass, is not sufficient to extend the range of validity of the effective theory beyond lo = me- 
One is therefore forced to use the complete calculation or fits to it. 

The outline of the paper is as follows. In section II, we give the essential steps of the 
calculation, collecting our hnal formulae in an appendix. Numerical results are presented 
in section III, together with the large rUe expansion and the hts. Finally, in section IV, we 
consider some of the possible astrophysical and cosmological implications of our results. 


II. COMPUTATION OF THE AMPLITUDE 

The leading Standard Model (SM) diagrams contributing to the processes in eq. 
are given in £g. 1, where permutations of the photon legs are understood. Momenta and 
polarization vectors are denoted by the sets {z}, {j} and {k} (z, j, fc = 1,2,3), with 
{1} = {pi, e“(pi,A)}, {2} = {p 2 , e^(P 2 ,p)} and {3} = {ps, e^(p 3 ,a)}, where A,p,a = ±. 
All momenta are incoming and q, flowing in the direction of the fermionic arrow, is the 
virtual integration momentum. 

{j} (a) 

; z 

P4 Z/ P5 

FIG. 1. SM leading diagrams contributing to five-leg photon-neutrino processes. 

At low energies, far from the scale of confinement, the leading contribution is given by 
diagrams involving only electrons in the fermionic loop. It is precisely the appearance of 
rUe as a scale, instead of M\y (which is the scale governing the four-leg photon-neutrino 
reactions), that makes such five-leg processes relevant already at energies of the order of a 
few MeV. 

By denoting with Aij^ and Bij^ the contributions coming from the diagrams (a) and ( 6 ) 
in fig. 1, the total amplitude M reads 

M(X, p,a) = [(Ai23 -|- A 321 ) -|- (Ai32 -|- A 231 ) -|- (A 213 -f A 312 )] 

+ [(-B123 + -B321) + (-B132 + -B231) + (-B213 + -B312)] , ( 2 ) 


(b) 

{j} 
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where, for example 


^i 23 = -E<r^/^"?Tr 

T=± 


^ 2 i = -E<rE^”9Tr 

T=± 


'I^W^(^)-^^{P2, p)-^^{pi, X)-^ 

^( 23 ) ^2 Vo V-l 

i^w^^ApiA)-^i{p2,p)-^i{P3,a) 

Vl Vo V -2 h;_( 23 ) 


( 3 ) 


In the previous equations and throughout the calculation, we used the notations 

Q±i = Q±i T me , = ^±(P) T me , 

Q±i Qo h: Pi 1 Q±{ij) Qo i Pi i Pj 5 f?0 Q 5 

-D±i = Q±i ■ Q±i 5 -0±(ij) = Q±(ij) ■ Q±(ij) ■ 


(4) 


Furthermore, = (1 ± 75 )/ 2 , and a dimensional regularization is used to compute each 
separately divergent diagram. Finally, the current F^ reads 

2 / 1 


rM = 




and 


Az = (P4 + Pbf -Ml, vf = Ve±ae, Ve = -]- + 2s^ , Og = ^ , 


(5) 


( 6 ) 


2 2 

where sw and cw are the sine and the cosine of the Weinberg angle, respectively. The reason 
why the terms in eq. are collected in pairs is because, when adding the two terms in each 
pair, and using the reversing invariance of the y-matrix traces, the 75 contribution cancels 
P]. For example 


^123 + ^321 — ~ 2 ner^ J dPq Tr 


7''7^/(P3, cr)^/(p 2 , p);^/(Pi, 

V(23) h;2 VO V-1 


(7) 


Exactly the same work can be applied to each pair of diagrams of type B. For instance, we 
get 


B 


123 


-B 321 — —2 Azc^F^ / d’^gTr 


7''7^/(P3, a)^/(p 2 , 

^{ 23 } ^2 Vo V -1 


Aw{q) ’ 

( 8 ) 


where A^iq) = {q + P 2 + P 3 + P^Y — M^^. Since we are interested in the large-(or 
equivalently low-energy) limit, we are allowed to substitute 0 

Az ~ —MI and A\[r ~ , (9) 
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so that the sum of the diagrams in eqs. (|^ and (|) gives 


^123 + ^321 + -Bi23 + -B 32 I — —2(ne + fill (Pl: P 2 , PS, ^1 Pi O') witll 



where the replacement ~ has to be performed also in the current 

Therefore, the total amplitude reads 


M(A,p,a) = - 2 (l + ne)r^ [If (pi, ps, P3, A, p, a) +/^(pi,p2, P3, A, p, a) 

+ h^(Pi,P2,P3,A,p,a)], (11) 


where and J 3 come from the remaining pairs of terms in eq. and can be obtained 
from Jf as follows 


P2^P3, ^{P2,p) ^ e{p3,a) 

P2^Pii e(p2,p) e(pi, A). 



( 12 ) 


The reduction of M(A, p, a) to scalar one-loop integrals can then be performed with the help 
of the technique described in ref. |^. The general philosophy of such a method is using the 7 
algebra in the traces to reconstruct the denominators appearing in the loop integrals, rather 
than making a more standard tensorial decomposition |P] . The algorithm can be iterated in 


such a way that only scalar and rank-one functions appear at the end of the reduction, at 


worst together with higher-rank two-point tensors. 

In this paper we content ourselves with giving the hrst step of this reduction, namely all 
results expressed only in terms of 

• scalar functions with 3 and 4 denominators, 

• rank-1 integrals with 3 and 4 denominators, 

• rank-2 integrals with 3 denominators, 

• rank-3 functions with 3 denominators. 

This already provides an important simplihcation with respect to the standard decomposi¬ 
tion, in that the computation of tensors such as 



(13) 


is completely avoided. 
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The key ingredient in the case at hand is a suitable choice of the polarization vectors 

m 

^{pi,X)=N + WxMsh] 

^{P2,P) = N + u!p:^2h:^i] 

/(P3, a) = N [|^i|^2|^3W<7 + Wa|^3|^2|^l] , (14) 

where 

^ = (4(ft.p.)(p(.p3)(p.-P3)) ’ = ^(1 ± ■ (15) 

Inserting eq. (0) in the expression for and using the identity 

hhi>iw^ + w^i>ihh = + w^i>ihh) + 2 (P2 ■ P3) [Qo - Qii], (16) 


gives 


^2(Pi,P2,P3, A,p,a) = -/f(pi,p3,p2, A,a,p) + J2 ^(pi,P2 ,P3, P, with 


<42^(Pi,P2,P3,P,c^) = 2(p2-^ 3)1^3 |y t^”?Tr 


b;(23) 




(17) 


With similar arguments we also get 


^3(Pi,P2,P3, A,p,a) 
<43 ^(Pi,P 2,P3, A,p) 


= -A^(P2,Pi,P 3,P, A,cr) + J3 (PiiP 2,P3, A,p) 


= 2 (pi-p 2 )iV" Md^gTr 


T 


Qi Q 


(13), 


1 1 
VO V-2 


with 


(18) 


Therefore, since and J!^ are directly expressed in terms of differences of three-point 
functions, If given in eq. (0) is the only master integral we need to decompose in terms 
of simpler tensorial structures. By inserting eq. ([T^ into eq. ([T0| ) for each combination of 
photon helicities, we obtain four different expressions. For example, when p = a = —A: 


A^(Pi,P2,P3, A, -A, -A) 


/ d^q 


Dq D-i D 2 D[2d,) 


kl^(A,-A,-A), 


(19) 


with 
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-A, -A) = ml ([/iQ(23)123231Qo231]A + [/iQ-il32Qol32321]A 
+ [/i231Qo23ig2321]A + [/i32ig223igo231]A) 

+ [/xg_ii 32 goi 32 g 2 i 23 g( 23 )]A + [/xg( 23 )i 23 g 2 i 32 goi 32 g_i]A, (20) 

and where we used the notation [/rij/c ■ ■ -Ja = ■ ■ • w'^]- 

To illustrate how the reduction works, take for example the last term in eq. (PU]). Since 
P 3 ■ P 3 = 0 ) we may rewrite 


[/ig(23)i23g2i32goi32g_i]A = -[/rg(23)i23ig232goi32g_i]A, ( 21 ) 

and, from the identity 

^2'^3'^2^0 = (-0(23) ~ D2)^2^0 + (Oq — D2);^3^o + {Dq + ml)^3^2 , (22) 

the hrst step of the denominator reconstruction immediately follows. As a second example, 
consider a term such as [/r ■ ■ ■ igol ■ ■ -(a- Since — 0_i), we immediately get 

[/i-■-igol ■ ■-(a = (Oo - 0_i) [/r-■-1 ■ ■-(a, etc. (23) 


The hnal result of such a procedure is given in the appendix. An important remark is in 
order here. To obtain compact expressions, we made a large use of the Kahane-Chisholm 
manipulations over 7 matrices |]M. Such identities are strictly four-dimensional, while we 


are, at the same time, using dimensional regularization. Our solution is splitting, before any 
trace manipulation, the n-dimensional integration momentum appearing in the traces as H 


q + q. 


(24) 


where q and q are the four-dimensional and e-dimensional components (e = n — 4), respec¬ 
tively, so that q-q = 0. The 7 algebra can then be safely performed in four dimensions, at the 
price of having additional terms. In fact, the splitting in eq. ( [2^) is equivalent to redehning 
ml —>■ ml — qf from the beginning. The net effect is then the appearance of extra integrals 
containing powers of q^ in the numerator, whenever mf. is present in the formulae reported 
in the appendix. The computation of such integrals in the limit e —0 is straightforward 
PJ. For example 


.TT^ 

Do D., D 2 D(23) ^ ^ ’ 


(26) 


A standard Passarino-Veltman decomposition of the simple remaining tensorial struc¬ 
tures in terms of scalar loop functions, concludes our calculation. We implemented the out- 
coming formulae in a Fortran code, performing the phase-space integration by Monte Carlo. 
Numerical results are reported in the next section. 
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Our formulae remain valid also when including all neutrino species. In this case, only 
the hrst diagram in £g. la contributes, at leading order in uj/ me, because the second one is 
suppressed by powers of uj/m^^T-. Therefore, the inclusion of all neutrinos can be achieved 
by simply replacing (1 + Ve) with (1 + 3ve) in eq. ©• However, we only considered Ug in 
our numerical results. 


III. RESULTS 


In this section, we give numerical results for the three reactions in eq. (|^) and present 
simple approximations to the full computation. Our main motivations are to assess the 
range of validity of the effective theory, and go beyond it. 

The total cross sections computed using the effective Euler-Heisenberg Lagrangian are 




a 


eff 


(77 


771 /) = 
7Z/I/) = 


^ 777 ) = 


262 

127575 

2144 

637875 

136 

91125 


\m, 


r — 


r — 


UJ 

me 


10 


10 




UJ 

m. 


10 


( 26 ) 


with k = Gi 7 ’me(l + Ue) Effective and exact computations are compared in 

figs. I, I and for the three cases. Furthermore, table I shows the ratio between exact 
cross section and for several values of uj/me- 

From the above hgures and numbers it is clear that, for all three processes, the effective 
theory is valid only when, roughly, oj/mg < 2 , as expected. At larger uj, the exact compu¬ 
tation predicts a softer energy dependence with respect to the (a;/me)^° behaviour given by 
the effective Lagrangian. 

All above results are in full agreement with those reported in ref. [^. Notice also that 
exact and effective predictions approach each other in the limit of vanishing uj, as it should 
be. This provides ns with an additional, strong numerical check on the correctness of our 
computation. 

Since the exact formulae are quite involved, we found it convenient to look for an ap¬ 
proximation above the point 10 /mg = 2. With this aim, we performed a large-rrie expansion 
of the integrals in eqs. ([I 0 |) and (|T^, to determine the OiX/rn^^ correction to the effective 
amplitude for 71 / —> 771 /. By iteratively applying the equation 


k'^ + 2q ■ k 


((q -h ky - ml) q^ - ml (^2 _ ml){{q + ky - ml) ’ 


(27) 
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and systematically discarding terms smaller than (9(l/mg), one is left with only contribntions 
of the kind 


ml 




(28) 


in which the external momentnm dependence in factorized ont from the integrals. In the 
above formnla, p* stands for any prodnct of s external momenta, and t, I, r, s are constrained 


by4 + r + f — 2Z + s = 0. Thns, when r + s = 0 the integral in eq. (pSf) is divergent, when 
r + s = 2 it is hnite and proportional to l/m^, and so on. 

By snmming all relevant contribntions, we fonnd 

+ Ve)‘^a^ / a; \ 3163 / u! 


ct(7z/ —>• 77^) = 


127575 TT^ 


262 


U! 

m. 


20 \m 


12N 


(29) 


The hrst term in the r.h.s. of eq. coincides with the resnlt given in eq. ([26|), so this is 
an extra check of onr compntation. The second term has the right sign bnt, as expected in 
ref. [ 0 , it does not give nnmerical predictions that are nsefnl to extend the effective theory 
beyond a; = mg. The only option is then to £t the cnrves in figs. ^ ^ and The resnlts of 
the hts are 


a{'yiy 771^) = > 77zz) x r 

X exp [2.13317 - 2.12629 log^(r) + 0.406718 log^(r) - 0.029852 log^(r)], 

cr(77 —>■ 7Z/P) = (t'^-^‘^( 77 —>■ 7i/p) X 

X exp [4.42122 + 0.343516 log2(r) - 0.114058 log^(r) + 0.0103219 log^(r)], 

Cr(z/P —> 777 ) = [yV —>■ 777 ) X 7-“6.57374 

X exp [5.27548 - 0.689808 log^(r) + 0.15014 log^(r) - 0.0123385 log^(r)], (30) 


where the effective cross sections are given in eq. (^ 6 ]), and r 
hts are valid in the energy range 1.7 < r < 100. 


uj/me- All the above 


IV. DISCUSSION AND CONCLUDING REMARKS 


The processes 1/7 —> 1 ^ 77 , —>■ 777 and 77 —are of potential interest in stellar 

evolntion and cosmology. The hrst two reactions can ahect the mean free path of nentrinos 
inside the snpernova core, while the last one is a possible cooling mechanism for hot objects 
. The relevance of snch processes depends on the size of the cross sections when varying 


the centre-of-mass energy a;, and may be dehnitively assessed only throngh complete and 
detailed Monte Carlo simnlations. Nevertheless, some simple considerations can be made, 
also in connection with speculations that recently appeared in the literature. 
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Basing their results on the assumption 
(t(7z/ ^ 771/) = (To 


ijj y 

1 Mevj 


(Jo = 10 cm^ 


(31) 


and on the data collected from supernova 1987A, the authors of ref. 0 htted the exponent 


7 in eq. (|M|) to be less than 8.4, for a; of the order of a few MeV. The physical requirement 
behind this is that neutrinos of a few MeV should immediately leave the supernova, so that 


their mean free path is constrained to be larger than 10 ^^ cm. 

The effective theory predicts 7 = 10, while, using the exact calculation, a softer energy 
dependence is observed in the region of interest (see £g. A fit to the exact curve gives 


7 ~ 3 for 1 MeV < a; < 10 MeV, thus confirming the expectations of ref. |12 


A second interesting quantity is the range of parameters for which the neutrino mean 
free path for such reactions is inside the supernova core, therefore affecting its dynamics. 


Always in ref. |T^ it was found, with the help of Monte Carlo simulations, that for several 
choices of temperature and chemical potential, and assuming the validity of the effective 
theory (7 = 10), this happens when ci; > 5 MeV. Since the exact results are now available, 
it would be of extreme interest to see how the above prediction is affected. More in general, 
we think that the reactions in eq. (|I|) should be included in supernova codes. 

In ref. 0, it was speculated that these processes could also have some relevance in 
cosmology. Consider, in fact, the mean number V of neutrino collisions, via the hrst reaction 
in eq. (|l|), in an expansion time t equal to the age of the Universe T^ : 


N = (t( 7 z/ ^ 'y'yi')n^ct , = neutrino number density 


(32) 


In the above formula, rij, and t can be rewritten in terms of the photon energy at thermal 
equilibrium (ci; ~ fcT), expressed in units of 10^° K. By denoting this quantity by Tio, we 
get 

= 1.6 X lO^^Tfo cm“^ , t = 2 s. (33) 

When N is large, the neutrinos are in thermal contact with matter and radiation, while, 
for V ~ 1 (namely a ~ 10“^^Th)^ cm^), the neutrinos decouple. By using the formula 
in eq. (p 6 D , the resulting decoupling temperature is Tio ~ 9.5, namely uj ~ 8.2 MeV, 
therefore outside the range of validity of the effective theory. By repeating the same analysis 
with the exact result, we found instead that N becomes of the order of 1 at a; ~ 1 GeV. 
At these energies, other processes enter the game; for instance reactions involving different 
leptons, quarks and hadrons inside the loop in £g. la. Therefore, the reaction 71/ —> 77Z/ 
is no longer the only relevant process. In conclusion, the hve-leg reactions in eq. ( 0 ) are 
unlikely to be important for a study of the neutrino decoupling temperature, contrary to 
what the effective theory seemed to suggest. 
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Appendix 


We list here the result of the decomposition of the master integral in eq. performed 

with the technique of ref. [^. The starting point is the formula 

A^(Pi,P 2,P3, A,p,a) = f (Tq ^ ^ -A^(A,p,ct), (34) 

J VqV_iV2JJ(23) 

with N given in eq. (plsl). 

By dehning {i ■ j) = Pi ■ Pj and m = nie, the four possible helicity configurations read 


• A = —p = —a: 

A^{X, -A, -A) = -4 (1 ■ 2)(2 ■ 3) Dq {[/iQol32Qol3]A + 2 (1 ■ 3)[/ig(23)132Q_i];,} 

+ (D(23) - D 2 ) {4m2 (1 ■ 2) [ [/r231Qo231]A - 2[pliQo2 - 2go)]A(l ■ 3)(2 ■ 3) ] 
-[/ig(23)123igo2312g_i]A - [/rgol32igo23igo2]A } 

+ {Do - D 2 ) {2 (1 ■ 3) [ 4 (1 ■ 2)(2 ■ 3)[/il(3go - go3)]A + [/i(32 - 23)gol231]A ] 

-2 [/ig(23)123gol32g_i]A(l ■ 3) + [/rgol32igo23gol3]A} 

-16m2 (1.2)(1 ■ 3)(2 ■ 3)(Do - D_^)[p2>2Q^]2^ 

+8 (1 ■ 2)(1 ■ 3)(2 ■ 3)(Do + m^) {[/il32]A - [pUS]^ + [/i321];, - [p231]x} 

+8m\l ■ 2)(1 ■ 3)(2 ■ 3) {-[pQol231]x + 2 (2 ■ 3)[pliQo2 - 2Qo)]x 
-2 (1 ■ 2)[/i231]A + [p2312Qo]x - 2 (1 ■ 3)[p32Q_,]x} 

-4m2(l ■3)(2-3)[/i32igo231]A. (35) 


\ = p = a: 


A, A) = 8m\l ■ 2)(1 ■ 3)(2 ■ 3){D^23) - D2)[p2{Qol - igo)]A 
-4 (1 ■ 3)(2 ■ 3)(Do - D 2 ) {[/igoig232ig(23)]A 
W [ [pl23Q2l]x - [pl321Qo\x + 2 (1 ■ 3)[p21Qo]x - 2 (1 ■ 2)[p3Qol]x ]} 

+4 (1 ■ 3)(2 ■ 3)(Do - ^- 1 ) {[/ig232igo2g_i]A 

W [ [p2312Qo]x + [p2Qol23]x + 2 (2 ■ 3) {[p21Q^23)]x - [pl21]x) - 2 (1 ■ 2)[pQo32]x ]} 
+16 (1 ■ 2)(1 ■ 3)(2 ■ 3){Do + m2)m2 {[ 3 (2 ■ 3) - 2 (1 ■ 2) ] [pl]x 
-2 {pl23]x - 2 [ (1 ■ 3) + (1 ■ 2) ] [p2]x} 

-8m\l ■ 2)(1 ■ 3)(2 ■ 3) {[p21Qo32]x - [pl23Qol]x + 2 (2 ■ 3){ [ - 4 (1 ■ 2) ] [pl]x 

+ [pQol2]x + [+gol3]A + [+12go]A - [+123]a} } + Am^l ■ 3)(2 ■ 3)[/r2igo3123]A . (36) 
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• X = p = —a: 

A^(A, A, -A) = -Wrn^l ■ 2)(1 ■ 3)(2 ■ 3){D^23) - D2)[pl2Qo]x 
-4 (1 ■ 3)(2 ■ 3)(Do - D 2 ) {[PQ 0 IQ 2 I 23 Q 2 U 
W [ [p321Q2l]x - [p21321]x - [p31321]x + 4(1 ■ 2)[pl3Qo]x ]} 

+4 (1 ■ 3)(2 ■ 3)(Do - D.,) {[/ig(23)123go2g-i]A 

W [ [p23Qol2]x + 2 (1 ■ 2){ [p23Q2]x - [p321]x - 4 [p23Qo]x + 2 (2 ■ 3)[pQo]x } ]} 
+16 (1 ■ 2)(1 ■ 3)(2 ■ 3)(Do + m^)7n^ {[pl23]x - (2 ■ 3)[/il]A} 

+8m\l ■ 2)(1 ■ 3)(2 ■ 3) {[p23Qol2]x - [+2312go]A + [+32igol]A - [+gol321]A 

+m^ [p{23 - 32)l]x + 2(1 ■ 3)[p32Q_,]x 

+2 (2 ■ 3) [ -[+gol2]A + 2 [pl2Qo]x - 2 (1 ■ 2)[p2]x ] 

+2 (1 ■ 2) [ -2 [p23Qo]x - [l^32Qo]x + 2 (2 ■ 3)[pl]x ]} 

+4m\l ■ 3)(2 ■ 3) {[/i3213go21]A - 2 (2 ■ 3)[p21Qol2]x} . (37) 

• A = —p = a: 

A'^iX, -A, A) = 4Do(l ■ 2)(2 ■ 3) {2 (1 ■ 3){[ + 2 (1 ■ 2) ] [p231]x 

+2 [ (1 ■ 2)[p32Qo]x + (2 ■ 3)[pQol2]x + (2 ■ 3)[/igol3]A ]} + [/igol32go31]A} 

+4 (D 2 - D(23))(1 ■ 2)(2 ■ 3) {4 (1 ■ 3)[ + (2 ■ 3) ] [pl2Qo]x - 2 (2 ■ 3)[pl2Qo31]x 

+2 (1 ■ 2){ 2 (2 ■ 3)[pl2Qo]x - 2 (2 ■ 3)[plQo2]x + [pl2Qo31]x + [+2132go]A } 
+[/iigo23i2go]A - [+goi32igo2]A + [/igoi32goi3]A} 

-8 (Do - D 2)(1 ■ 2)(2 ■ 3) {2 (1 ■ 3){[ + (2 ■ 3) ] [pl3Qo]x - (1 ■ 2)[p23Qo]x} 

+ (1 ■ 2)[ 4 (2 ■ 3)[pl3Qo]x + [+13go31]A ] - (2 ■ 3)[ [/xl3go21];, + [pl3Qo31]x ]} 

-2 (Do - D_i)(2 ■ 3) {4 (1 ■ 2)2[;i3go231]A + m2[/i213go231]A - m2[/i2312gol3]A 

+2 (2 ■ 3)[/igol3go213]A - 2 (1 ■ 2){4 (2 ■ 3)2[;xl2go]A + 4m^il ■ 3)[p32Qo]x 

+2 (2 ■ 3)[ 2 (1 ■ 3)[p32Qo]x - [+2312go]A - [+32igo2]A + [+3go231]A - [pQol3Qo2]x ] 

+ [/i3go2312go]A + [+go3igo231]A} + [/i2132igo231]A + [/i23igo2312go]A 
-[/.goi3go23i2go]A - [+go3igo2i32go]A} 
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-8 (Do + ■ 2)(2 ■ 3) {4 (1 ■ 2)(2 ■ 3){(1 ■ - [^i2]x) + [l^l32]x} 

+ (1 ■ 3) [ -4 (2 ■ + 2 (2 ■ 3)[fil32]x + m%2 [/il32]A - [f-i231]x) ] } 

+4 (2 ■ 3) {m\l ■ 3)[/i231Qo213]A + 2 (1 ■ 2f{4 (2 ■ 3)2([/il2Qo]A - [/^WIa) 

+2m\l ■ 3)([/x231]a + [^^32Qo]x) + 2 (2 ■ 3)([/il2Qo31]A + [/x2132Qo]a) + [/x231Qo231]a} 
+ (1 ■ 2){-4 (2 ■ 3)"[ [/il2Qo31]A - 2 (1 ■ 3)[fil2Qo]x ] 

+2 (1 ■ 3) [2m\2 ■ 3)(2 [/il2Qo]A + [^^Qol2]x + [^iQol3]x) 

W(m2[/i231]A - [/ilQo231]A + [^i231Qo2]x)] + m2([/i2312gol3]A + [/i23igo231]A)}} • 

(38) 


The expressions for I 2 and are given in terms of If in eqs. ( p!7D and ([T8|). 

Another advantage of the above redaction techniqne is the possibility of checking the 
calcnlation withont explicitly compnting the loop integrals. For example, with any arbitrary 
choice of Pi,P 2 ,P 3 ,(l and me = m, the expressions in eqs. ( 0 ) and (p5D shonld coincide 
nnmerically. 
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TABLES 


w/rrie 

—> 77ZV 

77 ^ ^vv 

vv 777 

0.3 

0.969(8) 

1.09(1) 

1 .20(1) 

0.4 

0.923(6) 

1.17(1) 

1.37(1) 

0.5 

0 .888(6) 

1.28(1) 

1 .68(1) 

0.6 

0.852(4) 

1.47(1) 

2 .20(1) 

0.7 

0.826(5) 

1.80(1) 

3.17(1) 

0.8 

0.811(5) 

2.41(1) 

5.31(2) 

0.9 

0.819(6) 

3.95(2) 

11.88(3) 

1.0 

0.880(7) 

23.1(1) 

176.3(3) 

1.1 

1.19(1) 

18.2(1) 

94.7(2) 

1.3 

1.71(2) 

8.31(5) 

31.6(1) 

1.5 

1.44(1) 

3.37(2) 

11.9(1) 

1.7 

0.996(8) 

1.40(1) 

4.96(3) 

1.9 

0.635(4) 

0.622(3) 

2.23(1) 

2.0 

0.503(3) 

0.424(2) 

1.54(1) 


TABLE 1. Ratio between exact and effective results for the three cross sections. 


The error on 


the last digit comes from the phase-space Monte Carlo integration. 
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